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The paper is devoted to the investigation of the accuracy of the 
Rozenberg approximation for the reflection function of a semi- 
infinite absorbing random media. We found that this approxi- 


mation can be used at y<0.5, where y= ./(1 — a)/(3(1 — g)). 


Here wp is the single scattering albedo and g is the asymmetry 
parameter. We also propose the modified approximation, which 
can be used at least up to y= 1-2, depending on the observation 
geometry. © 2002 Elsevier Science (USA) 
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INTRODUCTION 


The task of this paper is to study the reflection function of a 
semi-infinite turbid medium. This limit can be easily achieved 
experimentally. For this one should increase the geometrical 
thickness of a scattering layer or concentration of particles till the 
measured reflection function reaches its saturation level (1). 
Thus, the case considered here is of importance for the spec- 
troscopy of turbid media in general and for the spectroscopy of 
colloids and foams in particular. It is clear that even small ab- 
sorption in a single scattering event accumulates to a large and 
easily measured effect due to a random walk of photons in a 
semi-infinite scattering layer. Thus, the spectral measurements 
of the reflectance function provide us with the information on 
absorption coefficient of particles, which can be used, e.g., for 
the retrieval of their chemical composition. 

The solution of the spectroscopic problem is greatly simplified 
if analytical relations between measured reflection functions and 
absorption coefficients are used (2, 3). One of such relations has 
the following form (4-6), 


Rxo(E,n, p) = RYE, n, p)exp -U E, n, g)y), [1] 


where R..(&, 7, o) is the reflection function of a semi-infinite 
. . 0 . . 
absorbing turbid layer, R(E, n, g) is the same function 
for a nonabsorbing layer with the same microstructure, £ = 
cos vo, n = cos Ù, ¢ is the relative azimuthal angle between in- 
cident and observed beams, 7? is the observation angle, 1% is the 
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incidence angle, 


—4 1-9 [2] 
=Y 30-8)’ 


wo is the single scattering albedo, g is the asymmetry parameter, 


K K 
E40) a [3] 
1 
3 
Ko) = 5 f RÌ (E, nE +n) dn, [4] 
0 
1 Qn 
RE m) = 5 f R? E, n, o) dọ, [5] 
0 


The function Ko(€) can be approximately calculated with the 
following formula at € > 0.2 (7), 


3 
Ko(6)= 5 (1 + 26) [6] 


independently on a scattering medium microstructure. 

Equations [l]-[3] and [6] allow us to reduce the calculation 
of the reflection function of an absorbing medium to that at no 
absorption. This is an important result, which can be used for 
various spectroscopic applications (2). Equation [1] is called the 
Rozenberg equation. 

Our task is to study the accuracy of Eq. [1] for different ob- 
servation geometries and values of y. Also we slightly modify 
Eq. [1] to have a possibility to consider larger values of light 
absorption in a single scattering event. 


THE MODIFIED ROZENBERG EQUATION 


The reflection function is defined as the ratio of the intensity 
of light reflected from a medium in question to that of an ideally 
white Lambertian reflector (1). Clearly, this function for a semi- 
infinite turbid medium decreases with the parameter 6 = 1 — wo, 
which gives us the probability of photon absorption in a single 
interaction event. This follows from both Eq. [I] and Eq. [2] and 


0021-9797/02 $35.00 
© 2002 Elsevier Science (USA) 
All rights reserved. 


reflection function 


0 20 40 60 80 


incidence angle, degrees 


FIG. 1. 
dia with droplets characterized by the gamma particle size distribution 
f (a) = Aa® exp(—9a/a.y), where A=const, a.,,=6 um for different values 
of y = 0, 0.12, 0.37, 0.83, 1.17, and 2.56. It was assumed that the wave- 
length A is equal to 1.55 um. The refractive index of particles m = 1.3309 — 
ix, where x =0.0, 0.00001, 0.0001, 0.005, 0.001, 0.005. 


The reflection function Rə (£, 1) of a semi-infinite cloudy me- 


from the results of the exact solution of the radiative transfer 
equation, presented in Fig. 1 for different values of y. Naturally, 
it vanishes altogether if the probability of photon absorption, 
B is equal to one. In particular, we see that it is smaller than 
0.1 at the nadir observation and y > 2.56. The results presented 
in Fig. 1 are obtained solving the Ambartsumian’s nonlinear 
integral equation for the reflection function of a semi-infinite 
medium (8). 

It should be pointed out that reflection functions for turbid 
media with the same value of 6, but with different values of g 
differ. Larger values of g < 1 lead to larger y and, therefore, 
larger reduction of the reflection function Rə (see Eq. [1] and 
Fig. 1). This can be understood on the physical grounds as well. 
Indeed, the value of g tends to 1 for highly extended in the for- 
ward direction single scattering diagrams. This makes the escape 
of photons more difficult. Actually, both their travel distance and 
travel time increase before their escape from a turbid medium as 
g — 1. This, of course, will lead to the increase of light absorp- 
tion by a medium. Clearly, the reflection function will decrease 
in this case. 

It should be stressed that the reflection function Re. (E, Nn, o) 
of a semi-infinite nonabsorbing medium only weakly depends 
on its microstructure (9). Thus, reflection functions of absorb- 
ing semi-infinite turbid media with the same values of y (but, 
possibly, with different values of 6 and g) do not differ very 
much (see Eq. [1]). On the other hand, these turbid media could 
be quite different as far as their microstructure (size and shape 
of particles) is concerned. This simplifies the solution of the 
spectroscopic problem in a great extent. 
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To check the accuracy of Eq. [1], we plot the results of numer- 
ical calculations of the function 


EE „ R&C, 0.65, 0) 
~ U(1,0.65,0)  Ra(1, 0.65, 0) 


[7] 


against y in Fig. 2 for cloudy media of different microstruc- 
tures. The numerical code, developed by Mishchenko et al. (8), 
was used to find values of R? (1, 0.65, 0) and Rœ(1, 0.65, 0). 
The function U(1, 0.65, 0) was calculated with Eqs. [3] and [6]. 
Clearly, the relation [1] holds till y > 0.5. For larger y, however, 
symbols in Fig. 2 start to depart from the straight line F = y. 
This is especially pronounced at y > 1. Such a behavior is due to 
assumptions used in the derivation of Eq. [1], which is valid only 
for weakly absorbing media (4). Note that symbols for different 
sizes of particles at y > | still lay on the same line. This makes 
it possible to suggest that the function Ro(&, n, p) still depends 
primarily on the parameter y (and not separately on parameters 
wo and g or other characteristics of the phase function of the 
turbid medium under study). However, this dependence is more 
complicated than that presented in Eq. [1]. One obtains by fitting 
data in Fig. 2, 


F=y(1—cy), [8] 


where c = 0.05. 
The accuracy of this fitting procedure is quite high as it follows 
from Fig. 3. Thus, we have instead of Eq. [1]: 


RoE, n, 9) = RYE, n, g) exp(—y(1 — cy)UE, n, y). [9] 


This equation can be used to improve the accuracy of Eq. [1] 
at y > 0.5. Note that the convenient and accurate analytical 
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FIG. 2. The dependence F(y) (see Eq. [7]) at aş = 4, 6, and 16 um. Other 
input parameters are the same as in Fig. 1. 
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FIG. 3. The same as in Fig. 2 but for the dependence F(y(1 — cy)). 


expression for the reflection function R? (£, n, g) was given by 
Kokhanovsky (9). 

Equation [9] can be used to find the function y(A) from spec- 
tral measurements of the reflection function R..(A). Namely, 
we have 


1— /1-—4@c 
2c ? 


1. /R? 
=> hi- E 

U Roo 
Clearly, the function y(A) is related to the absorption coeffi- 
cient x (A) of substances inside small particles. For instance, we 


have for large compared to the wavelength weakly absorbing 
particles (6), 


F [10] 


where 


[11] 


P= Bx A)xef, [12] 
where xef = 27 aef/àÀ. The value of B does not depend on the 
absorption coefficient and size of particles. It can depend, how- 
ever, on their shape and the real part of the refractive index (6). 
Thus, we obtain with account for Eq. [2]: x (à) ~ yA). 

Note that the value of c was obtained by fitting exact ra- 
diative transfer calculations at a fixed observation geometry 
(E= 1, n=0.65). Performing the same fitting procedure for 
other values of 7 at € = 1, we found that the following simple 
approximate equation holds: 

c(n) = 0.05 — 0.197 + 0.2677. [13] 

Due to the reciprocity principle (1), the same equation is valid 
for the function c(é) at n = 1. The error of approximation, given 
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error, percent 
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FIG. 4. The error of Eq. [9] at y =0.12(x =0.00001), 0.83(x = 0.005), 
1.17(x = 0.001) as function of the incidence angle at nadir observation. Other 
input parameters are as in Fig. 1. 


by Eq. [9] with account for Eq. [13], is presented in Fig. 4 at 
various values of y. We see that the error is smaller than 3% for 
incidence angles smaller than 80° and nadir observation. Thus, 
Eq. [9] can be safely used at y < 1.17. For larger y, however, the 
error increases. This is illustrated in Fig. 5. We see that the error 
can reach —30% at y = 211. and o = 12°. Remarkably, there is 
an angular region (% = 50-60°), where the error of Eq. [9] is 
smaller than 7% even at y = 2.6 (see Fig. 5). This finding could 
be of importance for planning correspondent light scattering 
experiments. It follows from Fig. 5 that errors of Eq. [l] are 
larger than those of Eq. [9]. Thus, Eq. [9] should be preferably 
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FIG. 5. The same as in Fig. 4, but also for new (Eq. [9]) and old (Eq. [1]) 


approximations at y=0.12(x = 0.00001), 0.83(x = 0.005), 1.17(x = 0.001), 
y=2.11(x =0.001, a,= 16 um), y=2.6(x =0.005), as the function of the 
incidence angle at nadir observation. 
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used for the interpretation of experimental data. Note that Eq. [9] 
transforms to more simple Eq. [1] as y — 0. Thus, both formulae 
have comparable accuracies at small values of y. 


RELATIONSHIP BETWEEN THE DEGREE OF 
POLARIZATION AND THE DIFFUSE REFLECTANCE 


The aim of this section is to show how formulae presented 
above can be used for the interpretation of experimental results in 
light scattering media optics. For this, we take one example from 
a current research. Namely, Vitkin and Studinski (10) showed 
experimentally that the product of the degree of polarization in 
the backward direction p to the diffuse reflection coefficient r 
of selected human tissues of a different pigmentation is almost 
constant (p ~ 1/r with the correlation coefficient 0.72). Thus, 
we have 

pr=D, [14] 
where D ~ const. They have used the nadir illumination of 
skin samples, which are effectively semi-infinite light scattering 
media. 

Let us show how Eq. [14] can be obtained from results pre- 
sented above. For this we need to have the expression for the 
value of r, which is defined as follows (1): 


1 


27 
1 
rey=— | dp f ndn RG... 
0 


0 


[15] 


The approximate solution forr (€) can be obtained expanding the 
exponent in Eq. [9] with respect to y and performing analytical 
integration. The answer is 
r (£) = exp(—y(1 — cy) Ko(§)). [16] 

Note that the value of a (€) = 1 — r (&) gives us the portion of 
energy absorbed in the medium under its illumination along 
the angle arccos £. Taking the integral over incidence angle 
and following the same procedure as above, we obtain for 
the diffuse reflectance coefficient pọ under diffuse illumination: 
p = exp(—y(1 — cy)) approximately. The total absorption is 
proportional to 1 — p in this case. The ratio (1 — e)/(1 — wọ) 
can be interpreted as an average number of scattering events in 
the medium. 

Next we need the solution for the degree of polarization p. It 
was obtained by Kokhanovsky (11) as y > 0, 

PE, n, 9) = PHE, n, y) exp(yUE, n, p)), [17] 
where p? (E, n,o) = p(&, n, ¢, B=). This equation was de- 
rived, using the Rozenberg formula [1]. Better accuracy for 
larger y may be obtained if one applies Eq. [9]. Namely, us- 
ing the same approach as discussed by Kokhanovsky (11), we 
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FIG. 6. The degree of polarization as the function of the single scatter- 
ing albedo according to exact radiative transfer calculations, approximations 
[17] (old) and [18] (new). Input parameters were 0 = 6.7°, 3 = 32.5°, p = 0°, 
g=0.85, wo = 0.95, 0.97, 0.98, 0.99, 0.995, 0.998, 0.999 for Deirmendjian’s 
Cloud C1 phase function (12) at the wavelength 0.65 um. 


arrive to the following approximate solution: 


PE, n, 9) = PE, n, gexp( — cy)UE,n,g)). [18] 
We contrast results, derived from Eqs. [17] and [18] and the 
exact radiative transfer calculations in Fig. 6. We see that the 
use of Eq. [18] is preferable for larger values of y. 

Now, taking the product of r(1) and p(1, 1), we arrive at 
Eq. [14] with 


D = p? (1, 1) exp(S), [19] 
where 
E=ay(1 — cy) [20] 
and 
_ Ky) 


We used values of £ = 7 = 1 as in the experiment of Vitkin and 
Studinski (10) here. It follows that the value of D is not constant 
due to the variability of parameters p? (1, 1), œ, and y from one 
skin sample to another. This explains the scattering of exper- 
imental results around the average curve p(1, 1) =const/r(1). 
The scattering is rather small. Therefore, we conclude that & is 
a small number and 

D= p(l, 1). [22] 
Our estimations give us that the value of y varied in the range 
of 0.7-1.8 (except one point, where r ~ 0.1) in the experiment 
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of Vitkin and Studinski (10). We derived this result from their 
Fig. 4, using Eq. [9]. Also we estimated œ as close to 0.17 
for their experiment, where we used Eq. [6] at €=1 and 
the value R? (1, 1)= 1.134, given by Yanovitskii (13) for the 
Heney-Greenstein phase function at g = 0.9. This value of g 
is typical for various tissues. Thus, we have € ~ 0.1-0.3 and 
exp() © 1.1-1.3 with the average value of 1.2 over the whole 
interval of measurements. 
Finally, we note that the relation 


r) pC, 1)= p% (1, 1), [23] 


which is valid as y —> 0, can be easily checked experimentally. 


CONCLUSION 


In conclusion, we have studied here the accuracy of Eqs. [1] 
and [9] for the calculation of reflection functions of semi-infinite 
turbid media. It was found that Eq. [9] improves the accuracy 
of Eq. [1] at y > 0.5. We also found that the parameter c in 
Eq. [9] depends on the observation geometry. This dependence 
is given by Eq. [13] at nadir illumination of a scattering layer. 
Derived Eqs. [9] and [18] can be used for the interpretation 
of experiments in light scattering media optics as it was shown, 
applying them to experimental data of Vitkin and Studinski (11). 

The final remark is that the dependence p ~ 1/r is not a 
new result. It was etablished for the first time by Umow (13) 
and studied in detail by Ivanov (15, 16). Its application to the 
planetary optics is reviewed by Hapke (3). 

The decrease of polarization with the reflectance of the turbid 
medium can be easely explained on physical grounds. Indeed, 
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the high brightness of a turbid medium in reflected light is due 
to highly developed multiple light scattering, which randomizes 
both photon directions and polarization states. 
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